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We present magneto-transport calculations for YBa2Cus07_ ( 5 (YBCO) materials to show that the electron- 
like metallic chain state gives both the negative Hall effect and in-plane anisotropic large Nernst signal. We show 
that the inevitable presence of the metallic ID CuO chain layer lying between the Cu02 bilayers in YBCO ren- 
ders an electron-like Fermi surface in the doping range as wide as p = 0.05 to overdoping. With underdoping, 
a pseudogap opening in the Q1O2 state reduces its hole-carrier contribution, and therefore the net electron-like 
quasiparticles dominate the transport properties, and a negative Hall resistance commences. We also show that 
the observation of in-plane anisotropy in the Nernst signal - which was taken as a definite evidence of the elec- 
tronic 'nematic' pseudogap phase - is naturally explained by including the 'quasi-uniaxial' metallic chain state. 
Finally, we comment on how the chain state can also lead to electron-like quantum oscillations. 

PACS numbers: 71.10.Hf,71.18.+y,74.72.Kf,74.25.F- 



An understanding of the pairing mechanism in cuprate su- 
perconductors relies on the underlying complex Fermi sur- 
face (FS) properties, and how they are derived from the nor- 
mal state pseudogap (PG) phase. Important clues to the 
origin of the PG state have been put forward by recent 
Hall effect measurements which observed a negative Hall re- 
sistance below the PG temperature in YBCO systems. |[T3l 
This result received further supports from the negative See- 
beck coefficient, El [3 strongly enhanced Nernst effect. l4l[T4ll 
as well as identification of quantum oscillations RHTQl with 
cyclotron frequency indicative of electron-like FS. These 
stimulating discoveries, which according to the conven- 
tional Onsager-Lifshitz paradigm] 11] implies a closed elec- 
tron pocket, point toward the topological changes in the FS 
in the PG state when compared with its large metallic FS in 
the overdoped sample. |[T2l Such small electron pockets do not 
arise naturally from the band-structure calculations consider- 
ing the CuC>2 planes and their total area is inconsistent with 
the nominal doping concentration. fT3l 

To explain these results, candidate proposals for various 
symmetry-breaking patterns have been offered, leading to 
drastic FS reconstruction into multiple Fermi pockets of dif- 
ferent natures. l^ffl[^[T0l[T4l[T5l[T8l[T9l We would, however, 
expect to see signatures of these pockets in spectroscopies 
such as angle-resolved photoemission spectroscopy (ARPES) 
and scanning tunneling microscopy/spectroscopy (STM/S). 
Yet, extensive ARPES 01221 . and STM(23l studies on vari- 
ous cuprate materials consistently have so far been able to map 
out only one Fermi pocket or 'Fermi arc' - both would repre- 
sent the same FS segment if the intensity on the pocket's back 
side is suppressed by the coherence factors - which is again 
of hole-like character. (3j 0) Furthermore, to explain the large 
negative Hall coefficient, one requires to have electron-like 
quasiparticles dominating over the hole quasiparticles on the 
FS, which is prohibited in hole doped system due to Luttinger 
theorem. On the basis of this reasoning, the recent experi- 
mental evidence of a field induced charge ordering! 18] will 
also find difficulty in explaining the negative Hall coefficient 



even if it is assumed to induce electron pocket. [19 ] Taken to- 
gether, the electron-pocket scenarios of various density wave 
origins 01 031 033 are inadequate to explain the negative 
slope of the magneto-transport data. 

In this study, we introduce a fundamentally new perspec- 
tive - we show that the emergence of the electron-like FS 
in the PG state can be rationalized by realistically consider- 
ing the contributions of the much ignored metallic CuO chain 
bands in YBCO. The unambiguous existence of the chain 
states on the FS has been well established by many ARPES 
studies from extreme underdoped region (p = 0.05) to over- 
doping (p = 0.29 and above). [Q21 [23 Its contribution to c- 
axis transport l24l and quantum oscillation in specific heat] 25 ] 
measurements is also studied earlier, further supporting the 
metallic character of chain state. The chain states seem to 
exhibit neither the PG opening nor any considerable doping 
dependence. Nevertheless, it dominates in the electronic spec- 
tra in the PG state due to its interplay with the Cu02 states. 
Our reasoning that underlies this conclusion is that the PG 
opens in the Cu02 bilayers, truncating its full metallic FS into 
small segments. Irrespective of the specific origin and nature 
of the small FS, it contains much reduced hole-carrier con- 
centrations, determining the nominal hole doping of the sam- 
ple. Therefore, the low-energy electronic state becomes ef- 
fectively characterized by electron-like charge carriers of the 
chain state. Above the PG regime both in temperature (T) 
and doping (p), the large Cu02 state is recovered, restoring 
the predominant hole-carrier on the FS of YBCO. The re- 
sulting two-carriers FSs quantify the negative Hall effect, (T3J 
and also the enhanced Nernst signal with prominent in-plane 
anisotropy (H [14]| due to the 'quasi-uniaxial' character of the 
chain state. 

Tri-layer lattice model:- Encouraged by the ARPES 
results 021 [22| and the first-principle calculations. fT3l we for- 
mulate a unified non-interacting three band model, originating 
from 2D Cu02 bilayers and an ID CuO chain. (21 [27) To effec- 
tively model the reminiscent nodal FS segment of the Cu02 
bands in the PG state, we consider some form of density-wave 
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FIG. 1. (Color online) (a) Paramagnetic dispersion obtained from 
the eigenstates of the Hamiltonian in Eq.[T] Blue lines are the bond- 
ing and antibonding bands arising from the C11O2 planes, while ma- 
genta line depicts the chain state. The light blue and red shadings 
represent the hole-like and electron-like character of the bands, re- 
spectively, (b) Corresponding computed FSs are compared with the 
ARPES data[ 12 1 (plotted only in one quadrant of the Brillouin zone 
for visualization), (c) Same as in (a), but computed in the PG state 
in an underdoped region. The width of each line gives the associated 
spectral- weight, determined by the SDW coherence factors, (d) FSs 
in the PG state host only tiny hole pockets centering nodal points. 



order whose front side would also accommodate the 'Fermi 
arc' property. Constrained by this assumption, a Hartree-Fock 
like order parameter is introduced to the Hamiltonian with a 
periodic modulation Q = (ir, 7r), which is commensurate with 
the CuC>2 lattice plane. We can think of a residual spin- (3) or 
charge-density wave (18) (S/CDW) or a d-density wave (4) or 
some other commensurate FS instability as the microscopic 
origin of such ordering, however, the macroscopic properties 
we aim to describe here do not rely on these details. On 
the basis of these experimental considerations, we deduce the 
general form of the Hamiltonian in the Nambu representation 
per unit cell: 



h.c. 
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where c\ ka (cika) creates (annihilates) an electron with mo- 
mentum k and spin a on the i = p,p',c tri-layers. The 
singe-particle dispersions ^ are defined within tight-binding 
expansion of the Cu d^.^-orbitals (hybridized intrinsically 
with O p orbitals) hopping on a periodic lattice in the tetrago- 
nal crystal as given in Refl28l 

Next, we exemplify how a small hole-pocket forms in the 
PG state. To simplify the analysis, while retaining the salient 
spectroscopic features, we adopt a k -independent interaction 
potential V to be same for both bilayers states, while set- 
ting V = for the chain layer, as the latter do not exhibit 
any PG opening. We set the spin index a = —a to con- 
struct a SDW modulation. In this case, V encodes the on- 



site Hubbard [7, and a order parameter S, as Vi = USi, with 
&i = (4(k+Q)a c ika) d = P,P*\ where the canonical average 
is taken over the entire Brillouin zone. We use the value of ef- 
fective [7=1.59 eV, deduced from earlier calculations, (3| [32) 
and the order parameter S is evaluated self-consistent at each 
T and doping. See supplementary material 1 3 3] for details. We 
emphasize that our choice of a particular form of density wave 
does not preclude us from analyzing the full problem, rather 
improves the reliability of the results with respect to a single 
parameter [7. The choice of Q = (tt, tt) order however gives 
the FS topology which agrees well with most ARPES data. 

In Figs. [TJc) and (d), we notice that a direct band gap is 
opened at Ep everywhere in the Brillouin zone, except at the 
nodal points, giving rise to a hole-pocket in the CuC>2 states. 
The width of each line depicts the associated quasiparticle 
weight which is in detailed agreement with experiments, given 
that the weak intensity of the shadow bands can be detected 
by ARPES. l34l[35ll It is important to note that while the pa- 
rameters of the chain states are doping independent, its elec- 
tronic states posses a weak doping evolution via the interlayer 
coupling to the plane states. It is self-evident from the band 
curvature of the chain states at Ep that it continues to host 
electron-like FS (although open orbit) at all doping; see Fig. [2] 
for further illustration. The electron-like property of the chain 
state is also evident in the ARPES data of Ref.fHl 

Hall and Nernst effect calculation:- Relying on the afore- 
mentioned effective normal state Hamiltonian for YBCO, we 
now proceed to study the quasi-particle transport properties 
based on a semiclassical numerical calculation. Nernst exper- 
iment measures the transverse electric field, E, in response 
to a combination set-up of an externally-imposed temperature 
gradient, VT, and an orthogonal magnetic field, B.$4\ [14) 
The net electric current density J produced via this effect is 
essentially related to E and VT via linear response in-plane 
Hall conductivity tensor a and Nernst conductivity tensor a, 
as J=crE-aVT. In case of a perfect cancelation of the net 
charge current, we obtain a working formula for these two 
measured quantities: E = a~ 1 a\7T = — 0VT. Experimen- 
tally, the Hall resistance, p=cr~ l , and the Nernst coefficient, 
v=6/B, are measured independently via transport probes typ- 
ically by setting a weak magnetic field B=Bz along the c-axis 
of the lattice. 

In theory, a and a are calculated by solving the standard 
Boltzmann equations for the low-T DC transport. l9l [T0l After 
some tedious algebra we derive the formalism for these quan- 
tities in terms of the quasiparticle states and their associated 
weight (see supplementary material (33)): 



- ^E M K fc fi «Xfc sech2 liiEnk 



a 



k,n 



2 



k,n 



fiE n k 



(2) 



(3) 



Here E n k stands for the quasiparticle band with wavevector 
k, and band index n, of Hamiltonian in Eq. [T] Denoting the 
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corresponding Bloch eigenstate by ijj % nk , we express the trans- 
port matrix-element M^ k = ^nfe^nfc- M essentially projects 
the transport tensors from the Bloch quasiparticle representa- 
tion to the tri-layer lattice indices z, j. v^ k = dE nk /hdk^ is 
the quasiparticle velocity along the fi = x, y, z axes, and f3 = 
1 //c#T, where ks is Boltzmann constant. The differential op- 
erator that drives the system into a non-equilibrium state in re- 
sponse to the applied field is evaluated within relaxation-time 
approximation! 12]: Q nk = [-^(v nk x B) • V k +r k x ], 
where the relaxation time is taken to be constant both in k 
and T.[37 ] e, H and c have their usual meanings. 

The accurate evaluation of the T-dependence of the trans- 
port properties requires an experimentally constrained In- 
dependent gap parameters, in other words, an accurate de- 
scription of the T-evolution of the FS topology. It is well 
known that a mean-field order parameter overestimates the 
PG transition temperature. (3] |9) To overcome this difficulty, 
we adopt a phenomenological form of the T-dependence of 
the gap parameter V(T)=Vy/l - T/T S M\M with the onset 
T for the density wave order is set to be T s « 55 K, in close 
agreement with Neutron measurement for Nd-doped La-based 
cuprate|6|. However, at each T, the Fermi level is evaluated 
self-consistently to preserve the hole counting. 

Results:- Fig. [2] displays the T— evolution of the Hall co- 
efficient Rh = Pyx/B and Nernst coefficient v yx /T at two 
dopings at which the normal state Hall coefficient data are 
available (the value of B is scaled out in the presented re- 
sults). A visual comparison between theory and experiment 
reveals a systematic agreement between them. To relate the 
sign and magnitude of Rh and v yx /T to their corresponding 
charge carrier features, we introduce a function called 'quasi- 
charge weight', in analogy with the quasi-particle weight: 
Cij(k,u) = E F ) = En^^fcCfc^l^w - E nk ), where 
e nk = ±e at electron and hole like kp obtained from the slope 
of the bands at Ep, respectively, for n th band. We empha- 
size that in the orthorhombic real-space unit cell, the Hall cur- 
rent arises due to finite tunneling between two adjacent chain 
layers (which gives a warped chain FS in the /c-space). We 
plot the total the 'quasi-charge weight' in Fig. [2^ to d for sev- 
eral representative cases where the Hall coefficient is negative, 
zero and positive as a function of T and doping. 

The obtained results confirm our aforementioned postulates 
that at low-T when the strong PG shrinks the CuC>2 states 
into small hole-pockets, the dominant chain states produce a 
net electron-like quasiparticle character on the FS. With in- 
creasing T when the strength of the PG reduces, the area of 
the hole-pocket gradually increases, while the chain state re- 
mains very much unperturbed (note that with increasing hole- 
pocket area, the self-consistency allows the associated coher- 
ence factors to adjust themselves to maintain the same doping 
concentrations). At some intermediate temperature, T m , the 
'quasi-hole' and 'quasi-electron' weights exactly cancel each 
other to make Rh vanishes; see Fig.[TJb). At T > T m , the 
larger hole-pocket reverses the sign of the Hall-resistance, in 
Fig. [ije). Finally, above the SDW temperature, the metallic 
state is restored and the Hall coefficient becomes featureless. 
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FIG. 2. (Color online) (a)-(d) 2D color plots of the total 'quasi- 
charge- weight' dj (fe, Ef) (averaged over all layers) plotted at four 
representative cases. The red-white-blue color scheme gives the 
electron-like quasiparticle to no quasiparticle to hole-like quasiparti- 
cle states at a given k on the FS. (a) 'Quasi-charge- weight' calculated 
at an extreme low-T region for doping p — 0.10, where the Hall 
coefficient, Rh is negative and Nernst coefficient, v/T, is small, 
demonstrating the dominance of the electron-like quasiparticle state 
on the FS. (b) is chosen at an intermediate T where 'quasi-electron' 
and 'quasi-hole' weight compensate each other to lead to zero Hall 
coefficient and enhanced Nernst signal. At a higher T in (c) SDW 
disappears and the 'quasi-hole' weight dominates. Here Hall and 
Nernst coefficients flip sign. The plot in (d) is similar to that in (a), 
but at a higher doping p = 0.14. It should be noted here that, un- 
like in a strict ambipolar state, the maximum of the Nernst signal 
shifts from the zero of Hall coefficient, due to overlapping FS topol- 
ogy and matrix-element effects, (e) The computed T— evolution of 
the Hall coefficient, Rh, is compared with corresponding experi- 
mental data| 13 1 at two representative dopings of YBCO. The excel- 
lent agreement below the SDW temperature can be clearly marked, 
(f) Corresponding computed T-evolution of the Nernst coefficient, 
v a = v xy and z/& = v yx , at same T as in (e). The experimen- 
tal data for YBCO (shown in inset) is available at a slightly differ- 
ent doping. 1 14 1 Despite a characteristically good agreement between 
theory and experiments, we expect discrepancies due to the relax- 
ation time approximation as mentioned in the text. Results of v a and 
Vb reveal a large anisotropy, which is T-dependent. The source is 
this anisotropy is the presence of the 'quasi-uniaxial' chain state, ly- 
ing along the y-axis of the crystal. The anisotropy is strongest when 
the Nernst signal attains its maximum, in good agreement with ex- 
perimental data taken from Ref. [14] 



The same analysis applies to the higher doping result in that 
doping reduces the PG order, rendering a smaller hole-pocket, 
and hence less transport contribution. The results are in excel- 
lent agreement with the corresponding experimental data |[T3lL 
plotted with symbols in Fig.|2je) and (f). 

The discussion of Nernst effect proceeds similarly, al- 
though the fundamental mechanism for the enhanced Nernst 
signal is somewhat different. Experimentally, it is well es- 
tablished that the Nernst signal changes sign from negative to 
positive as well as it becomes abruptly enhanced below the PG 
temperature, but much above the superconducting transition 
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temperature, T c ,(4l Q31 see inset to Fig. (2jf). Despite some 
proposals that there exists vortex-type excitations even above 
T c which can drive Nernst signal] 7], the widely used ambipo- 
lar phenomena suggests that the compensating electron-like 
and hole-like FSs are primarily necessary to obtain an en- 
hanced 'normal state' Nernst signal. IHUlQl Within the present 
theory, such ambipolar state naturally evolves in the vicin- 
ity of the intermediate temperature T m , where 'quasi-hole' 
and 'quasi-electron weights' become comparable, and a pos- 
itive maximal of Nernst signal arises, as shown in Fig. [2jb). 
Quantum-oscillation measurement has also detected similar 
compensating electron- and hole-like FSs in YBCO.[8 ] Away 
from this T region, the positive Nernst signal diminishes 
as the underlying states depart from predominant ambipolar 
character, that is, when only one type of 'quasi-charge weight' 
dominates the excitation spectrum. Interestingly, even when 
the metallic state appears above the onset temperature of the 
density wave order, the so-called Sondheimer cancelation of 
the Nernst signal in a metallic state is still not strictly appli- 
cable here due to the coexistence of two-carrier charges on 
different bands. (42) In fact, the finite experimental value of 
the Nernst signal above T s , which remained an unexplained 
puzzle in the existing theories, is surprisingly well reproduced 
within the present theory. The results capture the salient ex- 
perimental features of YBCO data available at a slightly dif- 
ferent doping,| 14l[43) although some discrepancies are clearly 
visible. The adopted relaxation time approximation imposes 
limitation to quantitatively map out the detailed shape and 
peak positions of the Nernst signal. Furthermore, the present 
linear-response theory is mainly applicable to the low mag- 
netic field region; on the other hand, experimentally it requires 
sufficiently high field to expose the 'normal state' Nernst 
signal. Despite these approximations, the Hall and Nernst 
anomalies are robust features in YBCO, tied to the underly- 
ing electronic structure that is confirmed by ARPES lfT2l l22ll 
and first-principle results fT3l . 

In-plane anisotropy: -The added benefit of the present the- 
ory is that the existence of the the quasi- ID chain state natu- 
rally explains the associated in-plane anisotropy in the Nernst 
signal. For the ID state lying along the (lOO)-direction, the 
corresponding quasiparticle velocity is v x <C v y , which leads 
to a similar anisotropy in the Nernst conductivity a defined 
in Eq. |3j and hence in the Nernst coefficient v/T, plotted in 
Fig. [2ff). Most importantly, the in-plane anisotropy contin- 
ues to persist even above the PG phenomena disappears, both 
in experiment and theory, which adds confidence to our con- 
clusion that the PG phase in YBCO is unlikely to arise from 
any spontaneous rotational symmetry-breaking ordering. Fur- 
thermore, the metallic chain states have been demonstrated 
by ARPES data lfT2l to be present for doping range as large 
as p = 0.05 to p = 0.29 (and above), and all the exist- 
ing evidences of in-plane anisotropy via transport. I4ll441 neu- 
tron scattering measurement! 45 ] are obtained within this dop- 
ing range. This clearly indicates that the presence of chain 
state is responsible for the observed in-plane anisotropy in 
this system. It is worthwhile mentioning that most of the 



evidences for the electronic nematic phase in cuprate is ob- 
tained in YBCO family onry.l4l l44ll45l with some indication 
of it is presented in an analysis of the STM data in Ba-based 
compound. 1 46] 

Outlook and conclusion:- To further support our postu- 
lates, we refer to several other experimental facts where the 
present mechanism can offer a consistent explanation, at least 
in principle. (1) An important finding of the quantum oscil- 
lation measurements is that the cyclotron frequency of the 
election-like FS-which is proportional to the quasiparticle 
density at i^p-is very much doping independent in the op- 
timal doping region. O [TO) In this context, a phenomeno- 
logical model can be recalled which has shown that quan- 
tum oscillation can arise from open-orbit 'Fermi arc', by 
taking advantage of the electron-hole mixing of the Cooper 
pairs. 1 47] This argument finds support from an independent 
experimental consideration. [ 48 ] Based on this mechanism, the 
electron and hole quasiparticles, residing on different part of 
the FS in our present model, can also give rise to a quan- 
tum oscillation. If so, the doping independent cyclotron fre- 
quency arising from the weakly doping dependent chain state 
can naturally explain the experimental observation of the cy- 
clotron mass. Another explanation follows in that the open- 
orbit chain state becomes 'warped FS' and can create closed 
FS pockets if the orthorhombic-II unit cell is chosen, which 
within the conventional framework will give rise to quantum 
oscillation. (2) A second important result of recent quantum 
oscillations study is that the oscillation disappears at a crit- 
ical value of the doping for YBCO, and the cyclotron mass 
diverges as the critical value is approached from the high dop- 
ing side. Millis et fl/. fT5l have argued that the mass diver- 
gence is related to a Lifshitz transition where the multiple FS 
pockets connect to form an open (quasi- ID) FS. Our argu- 
ment follows similarly where we expect that in extreme un- 
derdoping, the Cu02 states become vanishingly small, and the 
residual ID chain states effectively cause a logarithmic diver- 
gence of the cyclotron mass. A rigorous calculation of the 
quantum oscillation would be of considerable benefit to con- 
firm our predictions. (3) Finally, we see no fundamental er- 
ror in postulating that the 'electron-like' and 'quasi-uniaxial' 
chain FS can also explain the negative Seebeck coefficient 
and thermopower signal, induced in-plane anisotropy in 
the spin- susceptibility measured by inelastic neutron scatter- 
ing study. (451 [49) Finally, we mention that since the chain 
state remains mettalic in double chain layered YBa2Cu40s+<5 
(Y124), |[T3l the above postulates continue to hold in this sys- 
tem, which is consistent with the observation of quantum os- 
cillation in Y 124. (50) 

Our realistic framework narrows the range of possible the- 
oretical models for the PG state in cuprates. It strongly sug- 
gests that the microscopic route toward understanding the na- 
ture of the PG physics may be generic to all cuprates con- 
taining a common Cu02 plane. Most of the other accom- 
panying normal state properties can then be explained as 
added features coming from the material specific crystallo- 
graphic differences, such as apical oxygens surrounding Cu 
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atoms on the C11O2 planes in La- and Bi-based supercon- 
ductors, chain layers in YBCO, chemically active Cl-ions in 
Ca2- iC Na ;E Cu02Cl2, among others. 

Note added: We argue that recent observations of charge 
density wave (CDW) do not alter our results. The reasons 
are: (1) The CDW has so far been observed in monolayer Bi- 
based cuprate l5Tl and YBCO[Q10 (52J, and its onset temper- 
ature and doping are very different from thePG temperature 
T* . On the other hand, the psuedogap phenomena and similar 
FS properties are obtained in all hole doped cuprates. (2) A 
theoretical study lfT9l has shows that a biaxial CDW induced 
nodal FS pocket is electron-like. According to the definition 
of electron-pocket, it is associated with a band bottom or band 
folding below E F , implying a gap opening along the nodal di- 
rection below Ep- But ARPES evidence for a gap along this 
line is yet far from definitive. Based on these facts, we argue 
that CDW is a secondary or weak material- specific feature, 
while the PG has more fundamental and generic origin. 

The work is supported by the U.S. DOE through the Office 
of Science (BES) and the LDRD Program and benefited by 
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SUPPLEMENTARY MATERIAL 



Eigenstates of the tri-layer Lattice model 



FIG. 3. Crystal structure and tight-binding hopping parameters of 
YBCU unit cell. 



Supplementary Fig. [3] presents the unit cell of YBCO crys- 
tal. The unit-cell contains two chain-layers and two Q1O2 
plane layers. In order to obtain a minimum low-energy model 
we assume that the inter-layer electron tunneling is only ac- 
tive between the the two nearest neighbor layers. In doing 
so, we are left with a tri-layer system which consists of two 
CuC>2 planes and one CuO chain, while the periodic bound- 
ary conditions on infinite lattice is imposed along all three 
dimensions. Furthermore, the chain state does not exhibit 
any gap opening through out the doping range of YBCO. The 
plane states undergo FS reconstruction which can be captured 
within a commensurate spin-density wave with modulation 
vector Q = (tt, 7t). The commensurate modulation doubles 
the unit cell in real-space. Using the standard Nambu notation, 
we define the tri-layer eigenfunction in the magnetic zone as 

^1 = c pfcf C p'fcf C cfcf C p(fc+QH' C p'(fc+QH' C c(fc+Q)| • 
In this notation the Hamiltonian presented in the main text in 
Eq. 1 becomes a 6x6 matrix H = ^2 k ^i^^: 



( £>ppk £>pp'k ^cpk Vppk ^ 

^pp'k ^ppk £,cpk Vp'p'k 

£,cpk £>cpk £>cck 

Vppk (,pp(k+Q) £pp'(fc+Q) (,cp(k+Q) 

Vp'p'k (,pp'{k+Q) £,pp(k+Q) (,cp{k+Q) 

(, C p{k+Q) (,cp{k+Q) (,cc{k+Q) ) 



(4) 



In obtaining the non-interacting dispersions we as- 
sume that Cu d x 2_ y 2 orbital contributes to the low-energy 
scales of present interest (as commonly used for all cuprates). 
The Cu d x 2_ y 2 electrons hop to their neighbors via O p x ^ y or- 
bitals which can be taken into account within effective tight- 
binding formalism. Q] |2) The hoping terms are depicted in 
supplementary Fig. [3] which gives, 

£ppk = -2t((/) x + <f>y) ~ ^t'(j> x (j>y - 2t"((j)2x + fay) ~ Hp, 
£,cck 2t c (J)y /i c , 

(,pp'k ~ —2tpp' (cj) x — <j>y)(j>zlA') (,cp = ~2t C p(j) z i4 Ll (5) 

where (j) a/Jj =cos (ak^) with /i=x,y,z. /i p?c are the chemical 
potential for the plane and chain states, respectively, which 
encode relative onsite energy differences and other crystal ef- 



fects between the two levels. t,t\ and t" are the 1st, 2nd and 
3rd nearest neighbor (NN) hopping on the 2D CuC>2 plane, 
t c is the NN hopping on the ID chain aligned along ?/-axis, 
and tp P ' (t cp ) is the plane-plane (chain-plane) hopping along 
z— direction. 

The eigenvalues for the upper 3x3 interblock matrix is ob- 
tained as 

Eck = e fc + Eok (6) 
EAk = e t ~ ^ofc (7) 
E B k = Xk- (8) 

Here subscripts 'A' and 'B' stand for the anti-bonding and 
bonding states created by the bi-CuC>2 layers, where 'C de- 
notes the chain state. e± = (xjfe ±£cck)/2,Xk = ippk±ipp>k, 
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and E ok = y((e^) 2 + 2^ pfc . Of course, the weight of 
the each layers are mixed between each bands which can 
be determined by their corresponding eigenvectors: u k = 
\ [e k ± E ok ] /icph- The eigenvalues of the full Hamiltonian 
can now be expressed as two split bands of the above three 
bands gapped by SDW order parameter Vij (3] 01 : 



E^u — ick =*= ^Ck 



E- 



Ck 

± 
Ak 



E Bk 



— ^Bk 



± 



(9) 
(10) 

(11) 



where ^f k = (E^k ± -EUfc) /2 and the self-consistent order 
parameters Vi are given below. The eigenvectors for the final 
eigenstates given above can be represented in terms of SDW 
coherence factors: 

1 



4c = ^(i+wVfe~ fc r 



v 



03, = i 



2 

a ik- 



(12) 



Here i 
in the chain state, o? Ck 



A,B,C. Due to the absence of any gap opening 
1 and f3(j k = at all momentum. 
With the definitions of the eigenstates and the coherence fac- 
tors, we construct the final unitary matrix U k that ultimately 
diagonalize the total Hamiltoanin given above in Eq. 



u k = 



U 



U k &Ak 



u k a A k u k a A k 
—a Bk ® B k 





U k pAk 
—pBk 





u k /3 Ak 

$Bk 



-OLBk 
&Bk 




$Bk 





~U k pAk 
-U^pAk 
-$Bk 



u k a Ak 

&Bk 



pBk \ 

-0Bk 






®Bk 





(13) 



The self-consistent order parameters are defined as V A = 
US A and Vb = USb (and Vc = 0), where U is the onsite 
Hubbard interaction which is kept to be same for both Q1O2 
planes. S A /b are me two order parameters, p is the dopin 
concentration which are evaluated by solving the following 
two self-consistent equations: 
1 

TV 



P : 



E 

i=l,2,cr 



Si 



(n it - mi) 



(14) 



27V ^ i 
k 

= ^E U ijk U^ jik f(E jk ), 



(k+Q)l Cik l 



(15) 



where / is the Fermi function and Ej are the eigenstates listed 



Eq. TOpT where indices i=A,B and j=l-6. is the num 



Oxygen content 



6,0 6^ 




300 



Q 200 -E 



0.1 0.2 0.3 




0.1 0.2 0.3 
Hole doping 



FIG. 4. (a) Schematic phase diagram of cuprate. (b) The doping 
dependence of self-consistently evaluated order parameter S (defined 
in the text). 



where the thermal average is taken over the whole Brillouin 
zone. 

Tight-binding parameters and self-consistent order pa- 
rameters 

We find that the tight-binding parameters that 
give a good description of the FS, in agreement 
with ARPES and first-principle band- structure arel2l: 
(t, t', t", t", t c , V j ^c P )=(0.38, 0.0684,0.095,0.25,-0.01,- 
0.0075) eV. We have kept /i c = -0.87 eV to be doping 
independent, while fi p is computed self-consistently. The 
band structure fitting is done at the overdoped sample of 
p = 0.29, shown in Fig. 1 in the main text. 

The calculation of self-consistency is carried out in the 
following steps. At a given doping p and onsite Hubbard 
(7=1.59 eV, we start our self-consistent cycle with an initial 
guess of chemical potential fi p and order parameters S A /b, 
and calculate doping p and S A /b by solving Eqs. (1) to 
(12) at T = 0. We iterate the entire loop until the self- 
consistent values of doping p, and order parameters S A /b 
converges. For the overdoped sample of p = 0.29, we get 
S A /b ~ an d /i p =-0.54 eV. At an extreme underdoped sam- 



, , —^1/25 ~ rp v — ■ ~ ■ • - — ~" 1 — 

ber of up-spin on the i m layer defined as ( c! (fc+Q)t c ifct ), p i es of p = 0.05 at which ARPES data are available for fitting 
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in Fig. 1, we get S A (T = 0) « 5 B (T = 0) = 0.3145 and 
li p = —0.4. At two dopings at which Fig. 2 is obtained in 
the main text, we find S A (T = 0) « S B (T = 0) = 0.25 
and /i p =-0.41 at T=0 to -0.43 eV at T = 55 K for p=0.10, 
and S A (T = 0) « S B (T = 0) = 0.22 and /i p =-0.38 at 7 =0 
to -0.41 eV at T = 55 K for p=0.14. The results are plotted 
in supplementary Fig.[4jb), and compared with the schematic 
phase diagram of YBCO in Fig. 1(a). Note that the temper- 
ature dependence of the order parameters is phenomenologi- 
cally assumed as Si = 5^(0) y/l - T/T s , where T s = 55 K 
is the SDW transition temperature. 01 



DETAILS OF HALL AND NERNST COEFFICIENT 
CALCULATIONS. 

The Nernst effect experiments measure the transverse elec- 
tric field response of a system to a combination set-up of 
an externally-imposed temperature gradient and an orthogo- 
nal magnetic field. The Nernst effect depends sensitively on 
anisotropies in the band structure. Similarly, the Hall effect is 
a powerful probe of the Fermi surface of a metal because of its 
sensitivity to the sign of charge carriers, which distinguishes 
between electrons and holes. These two probes have recently 
gained much attention in the context of cuprates, since they 
demonstrate very unusual behavior in the underdoped region 
when the pseudogap sets in. Apart from the superconducting 
fluctuation origin of enhanced Nernst effect in cuprates. 171 
purely quasiparticle picture is extensively proposed in var- 
ious contexts. [8HT0l The quasiparticle Nernst effect has 
been studied on the basis of the linearized Boltzmann equa- 
tion in the relaxation-time approximation. This is a reliable 
approach if the scattering rate is isotropic, since then the ne- 
glected 'scattering-in' contributions average out to zero. 

We start from the semiclassical Boltzmann transport equa- 
tion for quasiparticles (charge e) in a weak magnetic field B, 
driven out of equilibrium by a spatially uniform electric field 
E and temperature gradient VT. If we linearize the devia- 
tion of the fermion distribution function from its equilibrium 
state /(f fc ) in both E and VT as /'(&,) = / + 0fc> then the 
Boltzmann equation |[TT1l becomes: 

L tt° k g k = ^ Qk,k> (9k - 9k') (16) 

k' 



-eE 



Q k 



he 



(v k x B) • V fc . 



(17) 



(18) 



Here is the quasiparticle state whose velocity is v k = 
^ -1 Vfc£fc. The band index and the eigenvectors are implicit 
to the quasiparticle states in the above equations. Qk,k' is the 
scattering matrix-element for incoming state J2 k , Qk,k'9k' 
and outgoing state ^2 k , Qk,k'9k- For the present DC trans- 
port property, we assume the elastic scattering with Qk,k' = 
<K£fc — ik>)(lk,k>- Furthermore, we impose the standard 
relaxation-time approximation. lfT2l which implies that the 
system has enough time to loose its memory, i. e. scattering- 



in term is negligible. In that case, the scattering rate is gov- 
erned is entirely governed by the out-going state as ^ = 
9k Y^k' Qk,k', while relating its momentum dependence: 



1 

Tk 



^ Qk,k' = (qk,k') FS - 



(19) 



Within this approximation Eq. [T6|[T8] becomes 



9k = ^Vk 
V k = v k 
&k = — , 

he 



-eE - f, 



VT 



(v k x B) • Vfe ■ 



1 

Tk 



It is now convenient to introduce the band index n: 



ink 



d^nk 



(20) 
(21) 

(22) 
(23) 



From Eq. 23 the net electrical, J, and thermal current densi- 
ties, Q, that flow on the z £ Mayer (or orbital) can be calculated 
as from the diagonal term of the following matrix 



J* j = -2e^M%v nhgnh 
Q ij 



2 ^2 M n V nk£nk9nk, 
k.n 



where M is the matrix-element: 

Mii{k) = ^{kw n \k). 



(24) 
(25) 

(26) 



In Eq. 24 we project the densities from its band basis n (a 



band index is implicit on both sides of Eq. |2T ) to the layer 
indices, z, j. The factor '2' introduced in Eqs. 24||25 to count 
spin-degeneracy. 

Linear response theory: The interplay of electrical and 
thermal effects necessarily implies two conductivity tensors 
a, a, which relate charge current J to electric field, E and 
thermal gradient, VT vectors: 



a ZJ E-a lJ VT 



(27) 



The Nernst response is defined as the electricalfield induced 
by a thermal gradient in the absence of an electrical current, 
and is given in linear response by the relation E = —6VT. In 
absence of charge current (i.e. when J=0), Eq. (6) yields: 



E = (£ ij ') -1 a ij 'VT = -6 ij VT, 



(28) 



where Hall resistance, p=a -1 , and the Nernst coefficient, 
is=6/B, are measured independently via transport probes by 
setting a weak magnetic field B=Bz along the c-axis of the 
lattice. flUEHl In Eq. 



27 



both sides. From Eqs. 



21 



and 28 the layer indices z, j run on 



24 27 we easily deduce that 
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k.n 



dV 



nk 



k.n 



dE 

dVK 



nk 



<9(-VT) 



df(£nk) 
d^nk 

Of^nk) 
d^nk 



^- V M ij Q~ V serh 2 ^ ^ £ ' nfc 
o 2^ IVI nk V nk* L nk V nk secri 



k,n 

2~ 



k.n 



(29) 
(30) 



where we have 



r sech' 



used the identity 
In the usual manner. 



f g/(gnfe) A 

V K»h ) 



the differential 



operator Vt~ k can be arranged as a perturbative expansion 
in the magnetic field B$5\ [121 in order to obtain transport 
coefficients that do not depend on B. For this purpose we 
define Q nk = Y nk +A nk (Eq. [l8]), where T nk = l/r nk and 



A nk is the rest. Then 



ft; 



(9(£ 2 ). 



(31) 



The diagonal entries in Eq. 30 (in the /i, ^ basis, not in the z, j 



basis) are obtained from the zeroth order in B in Eq. 3 1 while 
the lowest-order contribution to the off-diagonal coefficients 
arises from the linear order in B in the expansion of Eq. 3 1 



To this accuracy, the expressions Eq. [30] can be simplified in 
form of the expressions 



a 



k,n 



/3e 3 B 
2hc 



k.n 



"nk dkll "nk Q k „ 



k.n 



f3 2 e 2 B 
2hc 



J2K j k<k^ 



k.n 



sech 



fiE nk 



v v OV nk jjl OV nk 



E nk sech 



2 / PE nk 



(32) 
(33) 
(34) 
(35) 



In a simplest model, we assume a momentum and tempera- 
ture independent value of relaxation time r k = To. A tempera- 
ture dependent value of To can improve the results, especially 
for the Nernst signal, however, with this simplest form, we 
are able to capture the salient features of both Hall coefficient 
and Nernst coefficient. The results in the main text are nor- 
malized to the constant value of B and To. The above expres- 
sions matches exactly with the one proposed earlier for vari- 
ous form of density wave order,[5, 9, 10] except the matrix- 
element terms. It is important to note that, if no other crys- 
tal symmetry is broken, then the Mott relation v xy = —v yx 
for Nernst coefficient can be recovered from the Boltzmann 
theory.[ 15 ] This implies that the choice of sign convention for 
the Nernst effect is arbitrary. To faciliate direct comparison, 
we have taken the same sign for | v a | and \v\\. 
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